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Abstract 

There is a well developed framework, the Black-Scholes theory, for the pricing of 
contracts based on the future prices of certain assets, called options. This theory 
assumes that the probability distribution of the returns of the underlying asset is a 
gaussian distribution. However, it is observed in the market that this hypothesis is 
flawed, leading to the introduction of a fudge factor, the so-called volatility smile. 
Therefore, it would be interesting to explore extensions of the Black-Scholes theory 
to non-gaussian distributions. In this contribution we provide an explicit formula 
for the price of an option when the distributions of the returns of the underlying 
asset is parametrized by an Edgeworth expansion, which allows for the introduction 
of higher independent moments of the probability distribution, namely skewness 
and kurtosis. We test our formula with options in the brazilian and american mar- 
kets, showing that the volatility smile can be reduced. We also check whether our 
approach leads to more efficient hedging strategies of these instruments. 



1 Introduction 



There are certain contracts called options which are negotiated in the stock 
market worldwide. An option is a contract that gives the buyer the right, but 
not the obligation, to buy or sell a given asset (a particular stock, an exchage 
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rate or even an index) at a future time T for a fixed price (or exchange rate 
or index), called the strike price. If the option is a buy option, its value Ct at 
expiration time T is given by: 

C T = Max(S T - K, 0) (1) 



where we denote S t the price of the underlying asset at time t and K is the 
strike price. The option price problem is to find a fair price C = C(S , t — 0) 
for the option today. 

These contracts are in principle intended to provide protection to the buyer 
against moves in the market that can result in large losses. In a seminal paper, 
Black and Scholes [1] showed that, if there is no possibility to make a profit 
without taking a certain amount of risk (the non-arbitrage hypothesis), then 
the price of contracts called options is given by the solution of a second order 
differential, the Black-Scholes equation. One of the basic assumptions of the 
model is that the returns of the underlying asset follow a geometric brownian 
motion: 

dS 

— = pidt + adW (2) 



where \x is the drift rate of the asset, a is called volatility and dW is a Wiener 
process with probability distribution of a normal distribution with zero average 
and variance a 2 = dt, that is N(0, dt). This stochastic differential equation can 
be solved resulting in: 

S T = 5oe (M-* 2 /2)T+aV™ (3) 



with W ~ iV(0, 1) so that, on average we have: 



( St ) = -— / dx S Q e^ 2 ^ T+ ^ x e~ x2 ' 2 = S e^ T (4) 



Given the process [2], the price of an option is given by the solution of the 
Black-Scholes equation: 

^+V^ + r^ rC-0 (5) 



with "initial" condition given by [1] and where r is the risk-free interest rate. 
The trained eye will recognize [5] as the backward Kolmogorov equation for 
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the conditional probability corresponding to the process [2] with \i substituted 
by r. 

An alternative way to find the price of an option is to use the Feynman-Kac 
theorem to write the solution of the Black-Scholes equation as an expectation 
value of the final price at expiration brought at present value by the risk-free 
interest rate: 

oo 

C = e~ rT J dxP(x,t)C T 

— oo 
oo 

= e~ rT J dx P(x, t)Ma,x(S T - K, 0) 

— oo 
oo 

= e~ rT J dx P(x, t)Max (S e^ 2 W+^x _ K , o) (6) 

— oo 

where the probability distribution in this case is symply a normal N(0, 1) dis- 
tribution. Solving the Black-Scholes differential equation [5] or the expectation 
value [6] yields the famous Black-Scholes pricing formula: 

C = S N(d 1 ) - Ke- rT N(d 2 ) (7) 



where N(x) is the commulative probability: 

,-z 2 /2 

7 2T 



r e~ z I 2 



and 



The basic assumption of the Black-Scholes framework is that the returns of the 
assets are normally distributed. However, this is not observed in the market. 
There is ample evidence that the real probability distributions of the returns 
present heavy tails and asymmetry. The market usually corrects for this by 
introducing a non-constant volatility, the so-called volatility smile. Figure 1 
shows schematically the relation between the volatility smile and the proba- 
bility distribution of the returns. If there is a higher probability of large moves 
in the tails, it is compensated by an increase in the volatility implying a larger 
price for the option. This volatility smile, however, bears no relation to the 
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Fig. 1. The non-gaussianity of the returns probability distribution is taken into 
account by introducing the so-called volatility smile. 



historical volatility obtained by a time series analysis and is simply a device 
to reproduce market data. 

In this contribution, we study an extension of the Black-Scholes framework 
for a non-gaussian distribution. There is a large number of non-gaussian prob- 
ability distributions that could be used. In particular, Lisa Borland recently 
obtained a solution for the option price in the context of the Tsallis distri- 
bution [2]. In a more general way, one could approximate the real probability 
distribution by an expansion around the normal distribution, as suggested by 
Jarrow and Rudd [3]. This is the so-called Edgeworth expansion. In section I 
we review the Edgeworth expansion, the distribution we chose to use to model 
the market returns. In section 2 we derive a closed form solution for the op- 
tion price using this distribution in a risk-free probability measure. In section 
3 we test our results with real market date and show that we can significantly 
reduce the volatility smile. We also perform a delta hedge test that shows no 
significant improvement over the traditional model. Section 5 concludes. 
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2 Edgeworth expansion 



In this section we will briefly review the derivation of the Edgeworth expan- 
sion. More details can be found in [4]. Let's assume a series of xi, x 2 , x n 
independent and identically distributed random variables with mean /i and 
finite variance a 2 . Defining the random variable 

X n = -f> (10) 
n 7^i 



the central limit theorem states that in the n — > oo, the random variable 
S n — y/n ^ Xn ~ IJ ^ approaches a normal distribution N(0, 1). However, we will 
be interested in the probability distribution before reaching this formal limit. 
This can be achieved by a cumulant expansion of the characteristic function 
of the distribution Xn(t) = E (e ltSn ^j, which for the normal distribution results 

in e~' 2 / 2 . The characteristic function can be expanded as: 



X n (t) = exp 



& 1 ,.N3 1 1 



(11) 



where we already used that E(S n ) — K\ — and Var(5' n ) = At 2 = 1. Performing 
a Taylor expansion in t of the characteristic function, collecting terms of the 
same order in n and performing the inverse Fourier transform in order to arrive 
at the probability distribution, we arrive at the Edgeworth expansion: 

( c K _ 3 a \ 

g {x) = 1 + ^{x 3 - 3x) + -^7-{x A - 6x 2 + 3) + i-(x & - 15a; 4 + 45x 2 - 15) z{x){!2) 



valid up to order 1/n and where we used the skewness £ = k 3 and kurtosis 
At = k.4 + 3, incorporating the factors of 1/n in these parameters. In the 
equation above z(x) is the gaussian distribution N(0, 1). 

The Edgeworth distribution is properly normalized due to properties of the 
Hermite polynomials but care must be taken in the parameters since it may 
not be positive definite or unimodal. In Figure 2 we show a graph of the 
allowed region. 



3 Option Pricing with the Edgeworth Expansion 



The distribution of returns of assets in the real markets are known to be non- 
gaussian, presenting heavy tails and asymmetry. We will model these returns 
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Fig. 2. Allowed range of the parameters for a positive definite and unimodal distri- 
bution. 

with the Edgeworth distribution. The advantage of using the Edgeworth dis- 
tribution is that it has a nice analytical form and we will see that it is possible 
to obtain a closed, albeit long, form for the fair price of an option on that 
asset. This facilitates the use of the results and the tests we will perform. 

The first thing we have to do is to find the risk-free measure in order to apply 
equation (6). In the risk- free measure the average price of the asset should 
obey: 

oo 

(S T ) = f dx g(x) So e^ 2 ^ T+ ^ x = S e rT (13) 

— oo 

This equation fixes the value of the parameter /i that must be used: 

^T = rT- log (l + ^(^v/T) 4 + ^-{aVff + f^v^) 6 ) (14) 

Notice that one recovers the usual result // = r when £ = and k = 3, that 
is, when we have a gaussian distribution. 

The option price is computed by evaluating: 

oo 

C Ed - = e~ rT J dx (?(x)Max (S e^ ! /^+-VT x _ ^ ^ (15) 

— oo 
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This integral can be performed yielding a closed form solution for the price of 
the option in the case of the Edgeworth distribution: 



C™»- = C$ a + (16) 

-S ({aVTfe + {v^T) A ex m + (<tVt) 3 (3(k - 3) + ?{x* m - 1) + 



gH-rT-xlt/2+aVTxr, 



72V2tt 

(aVf ) 2 (12£ - 3(« - 3)x rn + i 2 x m {x 2 m - 3)) + (aVf)(12^x m + 3(« - 3)(x 2 m - 1) + 
£V m - &£ + 3)))) + 

[f^^S - K) (3(« - 3)^(4, - 3) + l^a* - 1) + 



' e -rT-x^/2 



V 72V27T 

fx m (^-10< + 15))) + 
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^(^((aVT)^ - 3) + (av/T) 6 ^ + 12(orv^) 3 0) 



where x m = !2Ei^/jol_^ °" / 2 ) r j s the minimum value that the random variable 
can have for a non-zero integrand in equation (15). Notice that the usual Black- 
Scholes result is recovered when £ = and k = 3. In the next section we will 
use equation (17) in order to check if the volatility smile can be reduced. 



4 Edgeworth smile 



In this section we compare the Edgeworth model with the usual Black-Scholes 
type of model for the call options of the future dollar/real exchange rate 
negotiated at the brazilian BM&F and also call options of the S&P500 index 
negotiated at the Chicago Mercantile Exchange. The dollar/real option price 
are the daily opening prices in the market. 

We want to compare the volatitily smiles obtained from the Black-Scholes and 
the Edgeworth model. The volatility smile of the Black-Scholes model is cal- 
culated in the standard way, by numerically finding the volatility in the model 
that corresponds to the market price for a given a strike price. For the Edge- 
worth model, we first have to estimate the skewness and kurtosis parameter. 
This is achieved by minimizing the total quadratic error with respect to the 
model parameters (a, £, k): 

Mm J2 {PMarket ' PEd9e)2 (17) 

PEdge 



We make sure that the parameters are within the region of validity of the 
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Edgeworth expansion. Once £ and k are fixed by this procedure, we again 
numerically evaluate the volatility for different strikes in a given day. 

Some representative results are shown in figures 3. We can see that the smile 
is significantly reduced when using the Edgeworth distribution. 

In addition, we have also studied the performance of the Edgeworth model 
compared to the traditional Black-Scholes models in the delta hedge of port- 
folio, that is, if there is a real competitive advantage in using an alternative 
model in terms of making money. We concluded that the Edgeworth model 
does not present a quantitative gain in the delta hedge test as compared to 
the Black-Scholes model. 
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Fig. 3. Volatility smile for the Black-Scholes and Edgeworth distributions. These are 
smiles obtained from future Real/Dollar exchange rates with two different expiration 
dates. 
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5 Conclusion 



We have studied a non-gaussian extension of the traditional Black-Scholes 
model, using the Edgeworth expansion to parametrize the probability distri- 
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butions of the asset returns. We derived a closed form expression for the price 
of an european option on this asset. We verified with market data that this 
extension reduces significantly the volatility smile, implying that the model 
parameters are more robust than in the traditional approach. We also per- 
formed a delta hedge test of a given portfolio, where the extension did not 
show a quantitative improvement over the Black-Scholes model. This reduces 
somewhat the interest of large corporations in a practical use of this extension 
at this point. 

It would be very important to find a microscopic process that would lead to 
the Edgeworth distribution. This would provide a solid basis for performing 
Monte Carlo simulations for option pricing in this framework. 
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